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Abstract 

A system of three point vortices in an unbounded plane has a spe- 
cial family of self-similarly contracting or expanding solutions: during 
the motion, vortex triangle remains similar to the original one, while its 
area decreases (grows) at a constant rate. A contracting configuration 
brings three vortices to a single point in a finite time; this phenomenon 
known as vortex collapse is of principal importance for many-vortex sys- 
tems. Dynamics of close-to-collapse vortex configurations depends on the 
way the collapse conditions are violated. Using an effective potential rep- 
resentation, a detailed quantitative analysis of all the different types of 
near-collapse dynamics is performed when two of the vortices are identi- 
cal. We discuss time and length scales, emerging in the problem, and their 
behavior as the initial vortex triangle is approaching to an exact collapse 
configuration. Different types of critical behaviors, such as logarithmic 
or power-law divergences are exhibited, which emphasizes the importance 
of the way the collapse is approached. Period asymptotics for all singu- 
lar cases are presented as functions of the initial vortices configurations. 
Special features of passive particle mixing by a near-collapse fiows are 
illustrated numerically. 



1 Introduction 

The importance of point vortices in applications is due to the dominant role of 
coherent vortical structures in many 2D turbulent flows. Experimental and nu- 
merical studies performed in ||l|]-|§| have demonstrated, that in cases of driven or 
freely decaying 2D turbulence a number of concentrated vortices develops out of 
an originally unstructured flow. In many situations, dynamics of this flnite-size 
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vortices can be reasonably well approximated by point- vortex models ||9||-|[ll[, 
allowing the reduction of the study of a time dependent field, to a simpler Hamil- 
tonian system of N interacting particles. Point-vortex dynamics is an essential 
ingredient of punctuated Hamiltonian models, which have been successfully used 
to describe the evolution of 2D turbulence after the emergence of the vortices 
[|, 12, ISj. In these models, the advection of well-separated vortices is approxi- 
mated via Hamiltonian point-vortex dynamics; to account for the change in the 
vortex population toward smaller number of bigger vortices, dissipative merging 
processes (ij, |l^ are included for vortices which have approached each other 
closer than a certain critical distance. Point-vortex dynamics is responsible for 
bringing vortices together, i.e. it determines what kind of merger processes will 
occur and how often. 

The number of vortices in a fiow can be quite large, in which case a com- 
plete dynamical description gets intractable, but still can yield some important 
statistical quantities jl^, On the other hand, few- vortex systems can be in- 
vestigated in much more detail. Recent interest in this area was mainly directed 
towards the study of Lagrangian chaos in different settings ||l^-||2^. Apart from 
this, low-dimensional vortex dynamics is essential for the understanding of the 
evolution of many- vortex flows, since it describes their "elementary interactions" 
[|^ . On an unbounded plane one needs at least three vortices to get a non-trivial 
dynamics, leading to non-constant inter-vortex distances. On the other hand, 
motion of three vortices on an unbounded plane is always integrable, whereas a 
generic four- vortex system is already chaotic | p5| , p6| , [27| . Hence, a three- vortex 
system is the only generic nontrivial case, which allows a complete dynamical 
analysis. General classiflcation of different types of three vortex motion, as well 
as studies of special cases of particular interest were repeatedly addressed by 
many authors ||o|, |^ ||, ||, ||, |l||, see ||5| for a review. 

If we think about point vortices as an approximation of a finite-size patches of 
concentrated vorticity, than we have to be sure that they always stay separated 
by a distance several times larger than a characteristic size of the patch, in order 
for the approximation to be valid (see ||l^ for a direct comparison of finite- 
size patch dynamics with its point vortex counterpart). A natural question is 
how close can the vortices approach each other during their motion? A rather 
surprising answer is that three vortices can be brought to a single point in finite 
time by their mutual interaction. Aref jss) points out that this result was already 
known to Grobli more than a century ago. This phenomenon, known as point- 
vortex collapse was studied in j24|, |2^, ^ |3^ . The dynamics of exact collapse is 
relatively simple due to a self-similarity of vortex motion in this case: a triangle 
formed by vortices (vortex triangle) rotates and shrinks, but stays similar to 
the initial one. Vortex triangle area decreases linearly with time, at a rate 
determined by vortex strengths and triangle shape |Q . A spatial refiection of a 
vortex triangle is equivalent to a time reversal, so that a refiection of a collapse 
configuration will experience an infinite self-similar expansion with the same rate 
of change of triangle area. The finite-time collapse of three point vortices into 
one (or its reverse the splitting of one vortex into three) changes the topology of 
the fiow field, this feature itself has important consequences: this phenomenon 
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may play an important role in topologically driven transitions where vortices are 
involved, and it conceptually allows the concentration (collapse) or spreading 
(splitting) of vorticity even for non-viscous media. 

As a matter of fact, analysis of rare gas of vortex patches, performed in 
[ p8| indicates, that in the limit of low vortex density (vortex occupation less 
than one-two percent), strong interactions occur only between three vortices in 
a specific way, similar to the three point-vortex collapse. These processes may 
be thought of as resonance interactions in the vortex gas. 

For the collapse to happen certain conditions on vortex strengths and ini- 
tial positions have to be satisfied exactly, so that in a real system, probability 
of this event is zero: true resonances do not occur. However, actual vortices 
have some characteristic size, and when distance between them gets compara- 
ble with their size, vortices experience a considerable distortion, and may merge 
together, which means, that from a broader physical point of view, a collapse 
can be thought of as a process which brings vortices close enough (up to their 
size) to each other. For three initially well-separated vortices to collapse in the 
above sense, their strength and positions may satisfy "resonance conditions" 
only approximately. We refer to this kind of motion as near-collapse dynamics. 
Its importance stems from the fact, that inter- vortex distances are changed con- 
siderably (by orders of magnitude) during the motion, which brings the system 
to a different length-scale, where new physical mechanisms enter into play. 

The main goal of this article is to describe different scenarios of the 3- vortex 
system in near-collapse dynamics, in other words we perform a quantitative 
analysis of the motion in near-collapse situations, when the resonance initial 
conditions are sHghtly distorted, and/or the vortex strengths are not precisely 
tuned. These results are briefiy resumed in tables included in Section 3. 
We restrict our attention to the case when two of the vortices are identical, 
for which it is possible to map the dynamics of the vortices to a motion of a 
particle in a one dimensional potential. Our description is tailored to obtain the 
dynamical characteristics of the motion, such as the intervals of vortex distance 
variation, period of relative motion, etc, and even in a simplified situation of 
two equal-strength vortices, we have found a rich pattern of vortices dynamics. 
Shortly, depending on the parameters of the system, the vortices can perform 
periodic dynamics with different singular periods as well as periodic dynamic 
with non-singular period. More importantly we have found that the singularity 
of the periods can be of the logarithmic or pole type. This information about 
possible time scaHngs of the near collapse motion should be taken into account 
when many-vortex systems are considered and, particularly it is significant for 
studying advection in 3- vortex systems. We will come to this issue in a forth- 
coming publication. 

In Section 2 the basic equations of vortex motion, their symmetries and cor- 
responding first integrals are specified. Effective Hamiltonian for a three- vortex 
system with two identical vortices is introduced in Section 3. Using this Hamil- 
tonian, different types of near-collapse dynamics are classified and analyzed. 
In Section 4 we discuss various routes to collapse, i.e. the manner in which 
near-collapse motion types approach self-similarly contracting (expanding) so- 
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lutions as initial configuration is taken closer to the exact resonance. In Section 
5 we give an illustration of how the approach to collapse affects passive particle 
advection. 



2 Dynamical equations 

Point vortices are exact solutions of Euler equation (see for example which 
writes for the vorticity fl 

^ + [f^,V'] = o, (1) 

where [•, •] is the usual Poisson bracket, and ip is the stream function; for two 
dimensional motion of incompressible fluid fl = —V'^ip. Equation expresses 
the conservation of generaHzed vorticity along the path fines of the flow. 

A point vortex system is defined by a vorticity distribution given by a su- 
perposition of Dirac functions: 

N 

n{K,t)=J2k^S{x-x^{t)) , (2) 

1=1 

where x is a vector in the plane of the flow, ki is the circulation of i-th vortex, N 
is the total number of vortices, and Xi(t) is the vortex position at time t. Using 
this expression for the vorticity, and solving Poisson equation, in the Euler case, 
one obtains a stream function of a point vortex system. By Helmholtz theorem 
[|^, the motion of vortices is determined by the value of the velocity field at 
the position of the vortex. The equation of a point vortex motion is: 

at d-x.1 

where ez is the unit vector perpendicular to the fiow. The motion has an Hamil- 
tonian structure with the two Cartesian coordinates as conjugated variables, and 
by convention we shall drop i/' and use the notation H instead. The Hamiltonian 
H is given by: 

H^^Y.hk,U{x,,) (4) 

where Xij = |xi — Xj| and the interaction potential U(x) = — log(x) in case of 
an unbounded plane. The Hamiltonian (|4|) is invariant under translation and 
rotation, which implies both the conservation of the vortex momentum: 

N 

P = kiXi{t) — const. (5) 
1=1 

and vortex angular momentum Le^ with: 

N 

= fciX^(t) = const. (6) 
1=1 
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In this paper we consider a case of three vortices in an unbounded domain. 
We restrict our attention to the relative motion of vortices, taking inter- vortex 
distances Ri as prime variables. The notation we use is illustrated in Fig. |^. The 
invariance of the Hamiltonian under translations allows us a free choice of the 
coordinate origin, which we put to the center of vorticity (when it exists). Then, 
the other two constants of motion (energy and angular momentum) written in 
a frame independent form become: 

/ i/ = -^[fcifc2lni?3 + fcifc3lni?2 + fc3fc2lni?i] , . 

The equations of motion (||) yield the following non canonical system for the 
dynamics of the inter-vortex distances: 

k^^RiRi = AI-k(K^'^ - R^^) 

k2^R2R2 = A/Tr{Rf - R^^) (8) 
k^'R^Ra = A/7r(i?r' - i?2"') , 

where A is the area of the triangle A1A2A3 (see Fig. 1), and x refers to the 
time derivative of x, (See ^ for details). 

The physical system being defined, we shall now focus on the near collapse 
configuration. As mentioned earlier, for the three point vortex problem, under 
certain conditions, which depends both on the initial conditions, and the vortex 
strengths, the motion is self similar, leading either to the collapse of the three 
vortices in a finite time, or by time reversal, to an infinite expansion of the 
triangle formed by the vortices. These conditions of a collapse or an infinite 
expansion of the three point vortices are easily found. The first condition is 
immediate as for a collapse to occur K has to vanish. For the second condition 
we just need to look for conditions to obtain a scale invariant Hamiltonian. With 
such spirit let us divide all lengths Ri by a common factor A in the Hamiltonian 
H. We readily obtain H'{X) ^ H + (^ kikj) In A, we then obtain the collapse 
conditions: 

K = (9) 
Ef^O' (10) 

i. e, the harmonic mean of the vortex strengths ( p^ ) and the total angular mo- 
mentum in its frame free form (|9|), are both zero. Near a collapse configuration, 
the two conditions (^, ( p^ allow two different ways to approach the singular- 
ity. Namely, we can change initial conditions which changes the value of K , 
or change the vortex strength and modify the harmonic mean. Unfortunately 
the motion of the three vortices even though integrable is not easily computed 
analytically. It then seemed useful to restrict the problem to some less general 
case, assuming that the basic mechanisms of collapse should not differ much in 
the general case. To emphasize this last statement, let us first consider a general 
collapse situation, for which the total strength is ktot- The following conditions 
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are then fulfilled 

E^» = ^*°*'E^ = 0' (11) 

i i 

which are equivalent to 

^ fcj = hot ^ ktot^ . (12) 

i i 

The values for the vortex strengths corresponding to a collapse configuration 
are therefore located on the circle resulting from the intersection of the sphere 
of radius ktot^ and the plane imposing a total vortex strength equal to ktot- 
Since we have a circle, we shall assume that, from this point of view, no special 
configuration exists, and any should be sufficient to describe qualitatively the 
different possible behaviors. 

In the following we will be considering the case, when two of the three 
vortices are identical. Since the rescahng of time scales leads to a freedom in 
vortex strength normalization, we will assume 

fc2 = A:3 - 1 (13) 

Then, in order for a collapse to happen, the strength of the third vortex has to 
satisfy ( |To| ) which gives the resonance value of the first vortex: ki^ — —1/2. We 
will be mostly interested in the situations, when ki is negative and close to its 
resonance value, so we will denote 

k=\k,\, (14) 
and define the deviation from the resonance as 

5=l/2-k. (15) 



3 Near-collapse dynamics of vortices 

In the case when two out of three vortices are identical, there exists a convenient 
representation of the system dynamics in terms of a motion of a particle in a 
one-dimensional potential, parametrically depending on vortex initial condition. 
To derive this representation, we introduce a new set of variables: X = R\, 
Y = R2R3, and Z = R2 + R"^. The constants of motion (^ can be written as: 

A = e-^""" = Y''/X 

K = X-kZ. (16) 

where new parameter A is introduced instead of H in order to simplify formulas. 

An equation on X can be directly obtained from the first equation of (^, 
squaring it gives: 



^ - ^A^ — — — = ^A'k'^^ ■ (17) 
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Square of the vortex triangle area A can be found from geometrical identities 



A = Vl"|sin6l|/2 
X = Z -2VYcost 



(18) 



which leads to 

16A^ ^AY -{X- Z f . (19) 
Using the expression for the constants of motion ([l^), we obtain: 

^ -i^ ^2 ' (20) 

where Y = (AX)^/*^. This equation has a form of an energy conservation law for 
a particle of mass 1 and zero total energy moving in a potential V{X; A, K, k), 
defined by 

V{X,k,K,k) = g^2fc2y2 • (21) 



Indeed, equation (|2C|) can be rewritten as 

HeffiX,X;A,K,k) = P^/2 + V{X;AK,k) ^0, (22) 

with Hamiltonian equations 

X = dHef f /dP = P, P= -dHef / /dX . (23) 

Thus, dynamics of vortex configuration is governed by an effective Hamiltonian 
Heff, where the shape of the potential well V depends on the initial vortex 
positions through the values of first integrals A and K ( p^ , and the strength of 
the first vortex k. 

Effective Hamiltonian representation (p2|), gives a number of advantages, 
allowing to use simple standard techniques for one-dimensional conservative 
systems to find dynamical properties of vortex motion. In a sense, the problem 
is reduced to determining the shape of the potential V as a function of its 
parameters. 

Before proceeding to the analysis of the types of potentials, emerging in our 
problem, we have to mention the restrictions, imposed on the system by the fact 
that three vortices form a triangle. Since this was not taken into account in the 
derivation of the effective Hamiltonian ( p^ , it may put additional boundaries 
on the possible values of X. Note, that since the effective energy Heff in ( p2| ) 
is always zero, a motion exists only on the segments where the potential is 
negative. Therefore, triangle inequalities and the condition V < define the 
physical domain of X. Triangle inequalities 

\R2-R3\<Ri<R2 + R3, (24) 

written in terms of new variables can be obtained by taking the square of ( p^ : 
Z - 2VY < X < Z + 2Vy, which is equivalent to 

{X{k-l)+Kf-4:k^Y<0. (25) 
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The condition (|5|) is equivalent to the positiveness of the right hand side of 
equation ( p^ ) and translates that the square of the area of the triangle A1A2A3 
is positive. The condition V{X) < together with the triangle inequalities ( p5| ) 
imposes then 

{K - X f - Ak^Y > , (26) 

RecaUing that Z > 0, we have X > K, which together with inequality (p6|), 
gives finally 

X>K + 2kVY. (27) 

Now we have a precise information on where the motion should lie, and can 
proceed to the study of different regimes of near-collapse dynamics. Although 
the effective potential representation is valid for any A: 7^ 0, we will restrict our 
study (and our definition of "near-collapse") to the interval k e]0; 1[, since the 
special cases k = 1 (escaping vortex pair) and k = 2 (neutral tripole) , as well as 
the degenerate case k = (two vortices and a passive particle) bring up their 
specific singularities, unrelated to the collapse phenomenon. Let us consider 
different regimes such that each regime corresponds to a class of qualitatively 
similar motions [Q. Using the effective Hamiltonian representation (|2^), we 
infer that different regimes may be encountered corresponding to the number 
of roots of the potential V{X; A, K, k) lying within the physical domain. This 
number changes with variation of the parameters (A, K, k). 

Note, that the absolute value of K depends on the length units and can be 
scaled out of the problem. Indeed, the potential V and the physical domain 
boundary ( |27| ) are invariant under the scaHng transformation (compare to (p^) 

^^]ky A^A|if^ (28) 
and the problem can be reduced to the study of the three following cases 

K = l, K = -l, K = , (29) 

where the last case is singular. The different situations encountered have been 
condensed in the tables |[^, which are presented at the end of this Section. This 
is intended to briefiy summarize the result of this work and allow the reader to 
browse through this paper more easily. 

In the following we will discuss in details the different types of motions and 
critical situation encountered; most of the figures illustrating the problem will 
correspond to the values ( |29| ) , although we will keep K a.s a. scale parameter in 
formulas, keeping in mind, that only its sign is relevant to distinguish different 
regimes. We insist that the scaling of A in (p^), which can be rewritten in terms 
of vortex energy (^ as 

H ^ H + Slii\K\, (30) 

leaves A and H unchanged, when 6 = 0, i.e. when vortex strengths exactly 
satisfy the collapse condition ( p^ . The scales of the motion are determined by 
the value of the angular momentum K, while from a pure energetic point of 
view the motion would be thought as scale invariant. 
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3.1 Critical situations for 7^ 

To detect bifurcations, leading to the appearance of new roots of the potential, 
we look for the degenerate roots of (|2l|). As a first step, we just want to find 
these roots, leaving their detailed interpretation for the following paragraphs. 
Such roots Xc exist only for certain critical values Ac of the energy parameter 
A, and can be found from: 



(31) 



The above system can be easily solved if one notices, that the potential (21) is 
written in a factorized form: V = Vi* V2 / {Stt"^ k'^Y^) , with 

Vi = [{K -{l-k)Xf -APY], V2 = [{K - Xf - Ak^Y] (32) 

Then, if a solution Xc, Ac exists and is different from AT = 0, it should be either 
a double zero of Vi or V2 , or a zero of both Vi and V2 • 

We start from finding a double zero of V2 , which yields a system similar to 
where V is substituted by V2, that leads to: 

(X - K f - Ak^Y - , . 

2kX (X-K)- Ak^Y = ^ 

and after some substitutions we readily obtain, 

Xc, - K/25 

Ac, = 



(34) 



For Xci to be positive, K and 5 must have the same sign, so this bifurcation 
pertains to the physical region only for the cases K > k < 1/2 and K < 
0, k > 1/2 (it is easy to check that substitution of ( ^^ into (P^ ) turns the 
latter into an identity). Another solution, X — K, A = 0, corresponds to the 
situation when the negative vortex merges with one of the two others, reducing 
the system to a two-vortex case. Notice that X^ diverge as vortex strength 
collapse condition ( |lo| ) is approached {S ^ 0,k ^ 1/2), which announces a 
special behavior of the case 6 = 0, requiring a special treatment. 

Double zeros of Vi are found in the same manner. After some algebra we 
obtain: 

Xc2 = XcJ (1 - k) , , 

Ac, = (l + 2<5)Ac,/2 ^-^^^ 

and a solution X = K/{1 — fc), A = for the merged situation. Being propor- 
tional to Xci, Xc2 diverges in the same manner as (5 ^ 0, and requires sign of 
K and 6 to be the same, in order to He in a physical range. It also diverges 
when fc — > 1, as we have mentioned earher, this is the case of a scattering of a 
neutral vortex pair on a vortex, which is out of the scope of the present paper. 
Finally we note, that since Vi is proportional to the area of the vortex triangle 
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A, the critical equilibrium position X^^ corresponds to an aligned vortex config- 
uration; in a similar manner we deduce that corresponds to an equilateral 
configuration Q ^ . 

A third possibility, when both Vi = and V2 = 0, leads to the following 
system: 



This critical case also corresponds to an aligned configuration, where the nega- 
tive vortex of strength — fc is in the middle between the two identical vortices. 
In near collapse situation 2 — fc > 0, and since Xc^ has to be positive, this crit- 
ical situation appears only for positive total vortex angular momentum K > 0. 
Here a merged solution also exists {X = 0, A = +00), which corresponds to 
the merging of the two identical vortices. We also notice that the divergence of 
Xc3 for fc = 2 corresponds to a special case of a neutral tripole, when the total 
vorticity is zero, and the center of vorticity is not defined. 

Now we turn our attention to the three different specific cases distinguished 
by the sign of K. For each case, we will indicate which of the above found 
bifurcations do occur, and study the phase space trajectories of the effective 
Hamiltonian for different regimes. 



From the previous Section, we know that the number of critical energies for this 
case depends on the sign of S, i.e. whether fc>l/2orfc<l/2. A special case 
k = 1/2, 6 — will be treated separately. We start from the case k < 1/2, 
when all three critical cases (Q), (^5|), js^ ) belong to the physical region. If 
we keep the value of k fixed, and vary the energy, we will encounter various 
motion regimes, separated by the three critical energies. To illustrate the critical 
situations we plotted the potential as a function of X for the three different 
critical energies for fc = 0.2 in Fig. ^. The critical energy A^ corresponds to the 
maximum possible value for which motion can exist; for A = Ac2 an unstable 
equilibrium (saddle point) appears, and motion becomes aperiodic; and in the 
last case of A = Ac^ , another saddle point appears right on the border of the 
physical region (corresponding to unstable aligned configuration) . 

Motions regimes for different energy ranges are listed below (See Table |l|): 

1. A > Aci Motion impossible. 

2. Aci > A > Ac2 The potential has 2 zeros in the physical region, there 
exists a single type of periodic motion (see Fig. ||). When A Acj, the 
period diverges logarithmically 




(36) 



which yields: 




(37) 



3.2 The case is: > 0, A; ^ 1/2 



T(A) ^ l/21n|A- AeJ 



(38) 
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(see Fig. g) due to proximity of a saddle-point; motion acquires typical 
near-separatrix character of relatively short velocity pulses separated by 
long stays in the saddle-point vicinity. 

3. > A > The potential has 4 zeros and two different types of periodic 
motion exist (see Fig. . Their period diverges as 

T(A)^l/41n|A-AeJ, (39) 

when A approaches Ac2, Fig. |[ 

4. Ac3 > A The potential has 4 zeros and two different types of periodic 
motion exist (see Fig. |^). As A — ^ Acg (from both sides), the small-scale 
branch approaches a near-separatrix regime, its period diverges as 

r(A)^ln|A-A,3|, (40) 

for A = Ac3 this branch becomes aperiodic. 

For the case k > 1/2 only the critical value of X^.^ lies in the physical region, 
and only one regime of periodic motion exists (see Fig. |^) . 



3.3 The case < 0, /c ^ 1/2 

When the total vortex angular momentum is negative, the critical root ( ^ 
lies outside the physical region, and therefore we expect less variety of motion 
types. Otherwise, this case is similar to the K > case, in the sense that most 
phenomena described for K > Q occur, but in a reverse order; for instance, when 
fc > 1/2 we have three different regimes, which are analogue to the regimes of 
K > k < 1/2 case. The different regimes, illustrated by their phase portraits, 
are listed below (See Table 0): 

1. A < Aci Motion impossible. 

2. Ac2 < A < Aci The potential has 2 zeros and a periodic motion is possible 
(see Fig. ||). In vicinity of a saddle point (A « Acj) motion period diverges 
in the same way as in X > case. 

3. Ac2 < A The potential has 4 zeros and two different periodic motions are 
possible (see Fig. ||). 

On the other hand the situation when fc < 1/2, is analogous to the K > Q and 
fc > 1/2 situation. Namely, non of the critical values are physical, and we have 
only one type of periodic motion for all range of energies, see Fig. |l^. 

Now all the possible generic situations being described, we will proceed to 
the important singular cases X = or fc = 1/2. 
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3.4 The case K = 0. 



In this situation any rescaling of length does not affect the value of K, and it is 
only the value of A that controls motion scales in the system. The expression 
for the effective potential V{X) can be considerably simplified. It is convenient 
to introduce new variable 



then d21^) becomes 



ViU) ^ ^k'il-kr [U - ^—^ ] (U -I) . (42) 

The zeros are simple and the different conditions imposed on X imply a motion 
confined between 

^<U< ^ , (43) 

In this situation double zero do not occur. Note that for exact collapse case, 




k = 1/2, transformation to new variable (41) is singular and does not work. In 
this case the potential V{X) is constant; its value is negative only in a range 
of energies A S [1/2,1], which means, that collapse configurations have their 
energies confined to this interval. We will return to the discussion of collapse 
configurations in Section 4. 

In this singular case, to visualize the approach to collapse (A: ^ 1/2), it is 
better to have k fiuctuating and A fixed. This is due to the fact that, when 
k = 1/2, which corresponds to the merging conditions, the motion is possible 
only within a range of energies A2 > A > Ai. The phase portrait of trajectories 



taken for different values of k is shown in Fig. |11| and Fig. [IJ. These two figures 
are taken close to the values A = 1/2 and A = 1; we notice for the fc = 1/2 
case two straight lines corresponding to the two possible motions of expansion 
or collapse. Note, that when A ^ [1/2 , 1] a collapse cannot occur, which means 
that the trajectories for fc < 1/2 and fc > 1/2 are similar to those shown in 
Fig. |l^, but do not intersect on the top view, the results are summarized in 
Tableg. 

3.5 The case = 1/2 

This special case corresponds to "scale invariance", meaning that any length can 
be rescaled without changing the energy of the system. Therefore the energies 
of the critical situations should not depend on the value of K, which is scale 
dependent. The effective potential is now simply a fraction of polynomial and 
can be rewritten as 

V{X) = A (X - Xi) (X - X2) [X ~ X3) {X - Xi) /X^ , (44) 

where 

= ' = iJYTa ' = ' = TTa ' ^^^^ 



12 



and 

^ ^ (l/2 + A)(l/2-A)(l + A)(l-A) 

We have three critical values, the first two A^ = 1 and Ac2 — 1/2 correspond to 
the divergence of one root, while the third Ac^ = 1/4 corresponds to a double 
root situation X2 = Xr^. Depending on the sign of K we have the following 
allowed supports (See Table ^ : 

1. if > and A < A^g implies Xi > X > X2. The potential V{X) for 
the critical situation A = 1/4 is illustrated on Fig. |l^. As we approach 
collapse {K — > 0) this motion vanishes. 

2. K >0 and Acg < A < A^^ implies Xi > X > X3. The potential V{X) for 
the critical situation A = 1/2 is illustrated on Fig. |lj. As we approach 
collapse {K — > 0) this motion also vanishes. Nevertheless this case is 
interesting as an aperiodic motion very close to collapse is reached as 
A 1/2. The period of the motion as a function of |A — 1/2| is shown in 
Fig. |l^. It diverges according to a power law: 

(A- 1/2)"'/" , (47) 

(see the Appendix for its derivation). Phase space portraits of the trajec- 
tories are illustrated in Fig. |l^. 

3. A' > and Ac^ < A < A^, then X > X^. This motion is the closest 
to the collapse, as infinite expansion is possible, but X is bounded from 
below. In this situation we may see a collapse course, which rebounds on 
X3 and goes for an infinite expansion. As we approach collapse {K 0) 
X3 0, and depending on the sign of X a full collapse can occur. 

4. K > 0, Aci < A and K < 0, A < Ac^ motion impossible. 

5. AT < and Ac^ < A < A^ , then X > Xi. We have a similar situation as 
the one discussed for the K > Q and Ac^ < A < A^ . 

6. if < and A^ < A, then X-^ > X > Xi. Here we can anticipate a similar 
situation as the one discussed for the K > Q and Acg < A < A^^ . 

In this section all possible types of motion have been discussed, all the results 
for the different cases are summarized in the Tables |l[]^. A more detailed dis- 
cussion of the vortices behavior and trajectories, as the collapse configuration 
is approached, is presented in the next Section. 
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A 


iype oi motion 


Period 


K 


7 

k 


A > 


no motion 


no period 


K >0 


k < 1/2 


Aci > A > Ac, 


periodic motion 


T« l/21n|A- AcJ 


K >0 


k < 1/2 


> A > Ac3 


2 different periodic motions, 
1 large scale, 1 small scale 


1 l/4m|A — Acgj 


K > U 


k < 1/2 


Ac3 > A 


2 different periodic motions, 
1 large scale, 1 small scale 


T«ln|A-Ae3| 


K >0 


k < 1/2 


any A 


periodic motion 


no singularities for T 


K>0 


k > 1/2 


A < Aci 


no motion 


no period 


K <0 


k > 1/2 


Aci < A < Ac, 


periodic motion 


T« l/21n|A- Ac,| 


K <0 


k > 1/2 


Ac2 < A 


2 different periodic motions, 
1 large scale, 1 small scale 


l/41n|A- Ac,| 


K <0 


k > 1/2 


any A 


1 periodic motion 


no singularities for T 


K <0 


k < 1/2 



Table 1: Different types of motion for the general situation K ^ Q ,k ^ 1/2. 
The different critical values are 

= {25/Kf\ Ac, = i(l + 25){25IKf\ = i[(3 + 25)/Kf'. 



A 


Type of motion 


Period 


K 


k 


any A 


periodic motion, small scales 


period found numerically 


K = Q 


k < 1/2 


1/2 < A < 1 


collapse or expansion 


aperiodic motion 


K = 


k = 1/2 


any A 


periodic motion, large scales 


period found numerically 


K = 


k > 1/2 



Table 2: Different types of motion for the special situation K = ,k ^ 1/2. 
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A 


Type of motion 


Period 


K 


k 


A < Ac3 


periodic motion 


r«ln|A- Ae3| 


K>0 


k = 1/2 


Ac3 < A < Ac2 


periodic motion 


T« |A-AeJ-^/^ 


K>0 


k = 1/2 


Ac2 < A < Aci 


aperiodic motion 


no period * 


K>0 


k = 1/2 


Aci < A 


no motion 


no period 


K >Q 


k = 1/2 


A < A,, 


no motion 


no period 


K < 


k = 1/2 


Ac, < A < A,, 


aperiodic motion 


no period " '' 


K < 


k = 1/2 


Aci < A 


periodic motion 


T w |A- AeJ-3/2 


K <0 


k = 1/2 



Table 3: Different types of motion for the special situation k = 1/2 , K ^ 0. 
The different critical values are A^ = 1, Ac^ = 1/2, Ac3 = 1/4. 
* minimum approach X = X3 = K/ (1 — A). 
** minimum approach X = Xi = K/{l/2 - A). 
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4 Collapse scenarii 



As we have mentioned earlier, for the collapse to happen, two "resonance con- 
ditions" have to be satisfied. First, the sum of the vortex strengths inverses 
have to be zero. This condition does not impose any geometrical restrictions 
on vortex positions, for a given system of vortices it is either satisfied or not. 
Initial positions leading to collapse are specified by the second collapse condi- 
tion, K = 0. It also defines the range of energies, for which the self-similar 
dynamics can occur (for k = 1/2), but does not tell in which direction (collapse 
or expansion) it will go. 

For a given vortex configuration, let us consider a coordinate system with 
X-axis passing through the two positive vortices {ki — k2 — 1) and an origin in 
the middle between them, so that positive vortices have coordinates {—d/2,0) 
and (d/2,0), where d is the distance between them. Denoting the coordinates 
of the third vortex by {x, y) we can rewrite the condition K — Q as 

x'+y' = ^{l-k/2) (48) 

i.e. third vortex has to He on a "critical circle" centered at the midpoint between 
the two positive ones, with a radius (c?/2)-y/2/fc — 1. 

In the "resonant" case k = 1/2, the critical circle represents a set of initial 
conditions leading to a self-similar collapse or expansion. Level lines of A cross 
the circle in four points for A G (1/2, 1), see Figure |l^. These points are refiec- 
tions of each other in the coordinate axes, this symmetry lead the coordinate 
axis to divide the plane in four dynamically equivalent quadrants. A refiection 
of a vortex configuration is equivalent to changing the direction of time, so that 
points in adjacent quadrants have opposite directions of their dynamics. From 
the original equations (^), we deduce, that the parts of the critical circle lying 
in the I and III quadrants lead to a finite-time collapse, and those in II and IV 
quadrants to an infinite self-similar expansion. 

Four intersections of the critical circle with the coordinate axes are equilib- 
rium positions, where A reaches its limiting values (on the circle); the maximum 
value A = 1 corresponds to the equilateral triangle, and the minimum A = 1/2 
to collinear configuration. 

We define a rate of collapse as a rate of change of the squared distance be- 
tween the two positive vortices, i.e. as X. In this case, self-similar motion will 
have constant rate of collapse. Indeed, vortex velocities are inversely propor- 
tional to the inter- vortex distances, so that X — 2RiRi is independent of the 
scale of the motion. The rate of collapse is determined by the vortex energy: 

X{A) = v/4(l-A2)(l/4-A2), (49) 

it tends to zero when A approaches one of the equilibrium values A — 1/2, 1. 
Figure [T^ shows vortex trajectories for the fastest collapse case A = \/3/2. 

When the negative vortex is slightly away from the critical circle, motion 
type depends on whether A lies in the interval (1/2, 1) or not. If A G (1/2, 1), 
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vortices escape to infinity, their motion is aperiodic and unbounded. If the neg- 
ative vortex lies in the I or III quadrant, vortex triangle starts to contract in a 
nearly self-similar way, resembling an exact collapse trajectory, but by the time 
when X reaches its minimum value {K/{\/2 — A) for if < or K/{\ — A) for 
K > Q) the triangle evolves into a coUinear {K < 0) or isosceles {K > 0) config- 
uration, as shown in Figure |l^. After that vortices start to expand ad infinitum, 
asymptotically approaching a self-similar expanding trajectory. It follows, that 
a contracting self-similar motion in unstable, while an expanding one is asymp- 
totically stable, a result obtained in via linear analysis. When A ^ [1/2, 1] 
motion becomes periodic, X oscillates in an interval {K/ (1 — A); K/ (1/2 — A)) 
for if < or {K/il/2 - A); K/{1 - A)) for K > 0. 

In case k 1/2 unbounded motion does not exist anymore, yet the ratio of 
maximum vortex separation during the motion to the closest approach distance 
can be arbitrarily large as A: ^ 1/2. An interesting phenomenon is an existence 
of bounded aperiodic motions for particular values of vortex energy. They occur 
when an unstable equilibrium collinear configuration X = Xc2 appears at A = 
Ac2 . Apart from the logarithmic divergence of the motion period near Ac^ it also 
leads to a discontinuous dependence of the closest approach distance Xmin on 
initial vortex positions. Indeed, contracting vortices with A ~ — e where e in 
an arbitrarily small positive number, will not pass through X — Xc2, implying 
^min ~ Xc2 , while infinitcsimally different initial configuration with A = Aca + e 
will, leading to further contraction. This is illustrated in Figure |2^. 

To resume, approaching collapse means taking the limits K ^ and k 0.5. 
These Hmits are not commutative, and for certain situations, the result depends 
on how the limits are taken, e.g. these limits may induce the divergence of Xci 2 
or set them to zero or any value. Therefore a three vortex system can approach 
collapse in a number of different ways. For instance, if we take one limit after 
another, we can obtain the properties of the motion from the cases listed in 
Section 3.4 and Section 3.5. It is also interesting to study some of the scenarios, 
when the Hmits K 0^ and k — > 0.5+~ are taken simultaneously, leading 
to an interplay of different situations, described in the previous Section. As an 
example, consider a limit K ^ , fc — 0.5~ taken along the line K/S — C, 
where C is a constant. In this case, the saddle-point critical values {Xc2,Ac2 
tend to (C, 1/2), i.e. the scale of the inner potential well (bounded by ATcj) stay 
the same. At the same time the outer potential well spreads two infinity, and 
we have a situation, when two types of motion, corresponding to the same value 
of vortex energy, have entirely different length scales. 



5 Tracer advection near collapse 

As we have mentioned, a prominent feature of near collapse dynamics is a gener- 
ation of new length scales, which can differ from the length scale of the original 
configuration by orders of magnitude (see for instance (p7|)). This property has 
a strong infiuence on the mixing properties of the vortex fiow. In this Section 
we display some preHminary results on the difference of the advection of passive 
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particles (tracers) induced by collapse for a three-vortex system. We compare 
the motion of passive particles (tracers) in a velocity field of a near collapse 
configuration, with that in a typical "far from collapse" fiow. A stream function 
of the fiow, induced by moving vortices can be written as: 

1 ^ 

^Pi^,t)^-—Y,^n\^-Mt)\ (50) 

1=1 

corresponding equations of tracer motion are: 

drit) 



dt 



A Vi/) (51) 



where r(t) denotes tracer position vector. Vortex trajectories Xi(t) are periodic 
in a corotating reference frame, wher e (pl| ) have a structure of a Hamiltonian 
system with 1 ^ degrees of freedom j2l], |36[, so in a generic three- vortex fiow some 
of the tracer trajectories are regular, and some are chaotic. To visualize advec- 
tion patterns, we numerically construct Poincare section of tracer trajectories 
(in a corotating frame) , which is an expedient tool for estimation of the degree of 
chaotization of advection, location of the main structures in the chaotic regions, 
etc. In Fig. ^ an example of an advection pattern of a far from collapse fiow 
is presented. A large connected mixing region, regular cores around vortices, 
elliptic islands are typical elements of a tracer phase space in three- vortex fiows 
[3^ . Two most important length scales in the advection pattern, associated 
with its robust features, are the outer radius of the mixing region and the radii 
of vortex cores. Generally, these two scales are comparable (their ratio is about 
one order of magnitude or less) . However in advection patterns of near collapse 
vortex configurations, these length scales become considerably different, as the 
core radii tend to zero when approaching collapse see Fig. The figures Fig. ^ 
and Fig. are showed to emphasize that tracers chaotic dynamics in the near 
collapse case, is different than the dynamics without collapse. A more detailed 
analysis will be described in another article. 



6 Conclusion 

In this paper the behavior of a three point vortex system in a near collapse 
state was studied. We restricted ourselves to a special case where to vortex 
are identical which allowed the reduction of the system to an effective one di- 
mensional Hamiltonian study. We were then able to classify all typical types 
of motion which could occur in a near collapse case, which qualitatively agrees 
with previous classification of the motion of three vortices or collapse 

perturbations ||2^. On the other hand our classification also focus on dynam- 
ical properties of the system, which are exhibited through the representation 
of trajectories in the effective phase space, giving us for example a direct in- 
sight on the period of the motion for the closed trajectories. It is in this sense 
complementary of the previous general classifications. 
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Furthermore, the singularity of the vortex collapse is special in the sense 
that near collapse dynamics shows a number of nontrivial dynamical features. 
In particular, the possibility of vorticity concentration allows one to speculate on 
the fact that the collapse situation exhibits some basics mechanisms which are 
of importance in 2D turbulence and its models. For instance, it becomes clear, 
that three- vortex collisions play an important role in the aggregation of vorticity 
in 2D turbulent flows, especially at the late stages of freely decaying turbulence, 
when the vorticity occupation is rather low [Q. Among these three- vortex col- 
Hsions, near collapse ones are the more effective in bringing vortex together, and 
are thus the most important ones; this speculation is supported by the numer- 
ical simulations |Q. From this perspective, we give in this paper through the 
analysis of the effective potential, different typical length scales {Xc) and time 
scales (critical behaviors of the period), which shall be of importance in models 
attempting to describe 2D turbulence. For instance, the intricate structure of 
the vortex collapse singularity results in a very sensitive, non trivial dependence 
of the closest approach in the vicinity of collapse; it turns out that this is a 
discontinuous function of initial vortex positions. This is illustrated in Section 
4, by the flgure ^; in the same Section the conflguration of the three vortices in 
their closest approach has also been discussed, and a clear distinction is made 
between the coUinear and isoscele triangle configuration. These results could be 
of importance in deriving three- vortex merging rules analogous to the two- vortex 
ones used by current punctuated Hamiltonian models of 2D turbulence. 

From the dynamical point of view, the study of vortex motion near the pecu- 
Har critical situation of collapse, has also revealed itself interesting in exhibiting 
a variety of behaviors. Namely, we observed for some situations different peri- 
odic motions for which critical behavior of the growth of the period with respect 
to the control parameter A can be logarithmic, but it also can be as well a non- 
generic power-law growth with exponent 3/2 (See Section 3.5 or Table This 
last behavior is particular as in most cases, the approach of a singular critical 
point leads to a typical logarithmic divergence of the period, as for the near 
scattering situation discussed in |Q. 

Another motivation was how, the near-collapse dynamics influence the prop- 
erties of passive particles. Advection in the chaotic region of the flow due to 
three identical vortices is superdiffusive |^^; tracers statistics is non-Gaussian, 
and is characterized by power law tails in the tracer probability density function, 
and algebraic decay of Poincare recurrences statistics; both are the consequences 
of the presence of long regular flights in a chaotic tracer trajectory. We expect 
the advection in a general three vortex system to be anomalous too, but the 
question whether or not the power law exponents and self similar behavior of 
the the tracer distribution observed in [|9| are universal, remains open at the 
moment. With this in mind, it seems that the "extreme" dynamics encountered 
in near collapse dynamics, seem to be the most appropriate candidate to check 
if such universality exists. In this paper we limit ourselves to a short illustration 
of the advection patterns evolve as collapse is approached. We have observed 
, that as the system approaches collapse, the degree of chaotization increases: 
the chaotic region grows and the area of the regular island inside it shrinks. A 
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more detailed study of the dynamical and statistical properties of the advection 
in near-collapse flows is currently under way, and will be presented elsewhere. 

To conclude, we emphasize once again the fact, that near collapse dynamics 
comprises several types of motion. Their time and length scales differ consider- 
ably, depending on the way the collapse conditions are violated. 



Appendix: Power-law divergence 



From the effective Hamiltonian ( |22| ) the period of the motion in the case fc = 1/2 
and for instance K > and A < 1/2 (see Fig. |l^) is then simply deflned by 

T = 2 / , , (52) 

where X-y and X-^ are defined in (|4^). To obtain the asymptotic expression of 
the period as A ^ 1/2^ (where the — sign refers that we are taking a right limit 
A < 0.5), we use some simple equivalents, and write 

using the expression ( ^5|) for X\ and (ji^ ) for A, we obtain the asymptotic 
behavior of the period as a function of A and K 

3/1 



The period diverges as (1/2 — K)^"^!"^ when A approaches the collapse value of 
1/2, and if 7^ 0. Note that the period is Hnear in K. 
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Figure 1: Notation chosen for the probL 
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Figuro 2: The potential V{X) for the three critical cases, k = 0.2. For A = 
the motion appears, for A = two motions are possible, and for A = A^g an 
unstable aligned configuration is approached as X ^ X+ . 
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Figure 3: Phase space portrait close to for different energies for the case 
K = 1 and k = 0.2. The bottom figure is the top view of the upper one. Once 
the first critical situation is reached a motion is possible. 
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Figure 4: Phase space portrait close to for different energies for the case 
K = 1 and k = 0.2. The bottom figure is the top view of the upper one. Once 
the second critical situation is reached, the single possible trajectory reaches the 
separatrix, and the splitting in two possible motions occurs. 
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Figure 5: Period of the motion versus energy {AH = \H{X) — -ffc2l) for k = 0.2 
close to A = . We notice the logarithmic behavior of the divergence, we note 
also that the sum of the slopes of the two right branches (lower) is equal to the 
slope of the left branch (upper) . 
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Figure 6: Phase space portrait for A = Ac., for the case K = 1 and k = 0.2. A 
splitting occurs like for A = h.^^ but the triangle conditions prevent the existence 
of a third branch. We notice the sharp angle aX X = Xc^. 
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Figure 7: Phase space portrait for different values of A in the case K = 1 and 
k = 0.55. The bottom figure is the top view of the upper one. The motion 
remains periodic. As A increases, we notice the growth of the range of length 
scales explored by the motion. 
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Figure 8: Phase space portrait close to for different values of A in the case 
K = —1 and k = 0.55. The bottom figure is the top view of the upper one. 
Once the first critical situation is reached a motion is possible. 
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Figure 9: Phase space portrait close to for different values of A in the case 
K = —1 and k = 0.55. The bottom figure is the top view of the upper one. Once 
the second critical situation is reached, the single possible trajectory reaches the 
separatrix, and the splitting in two possible motions occurs. 
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Figure 10: Phase space portrait for different values of A in the case K = —1 
and k = 0.2. The bottom figure is the top view of the upper one. The value of 
A determines the typical length scale of the periodic motion (the value of K is 
fixed) . 
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Figure 11: Phase space portrait close to for different values of k in the case 
K = 0. The bottom figure is the top view of the upper one. For this singular 
case we notice the two lines at fc = 1/2 corresponding one to collapse and one 
to infinite expansion. On the projection plot (bottom), the trajectories k < 1/2 
and A; > 1/2 seem to intersect on those two fines. For A < cofiapse is 
not possible, and the trajectories do not intersect. We notice that close to the 
collapse condition k = 1/2, the values k < 1/2 describe a collapse motion, while 
fc > 1/2 describes the expansion motion. 
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Figure 12: Phase space portrait close to for different values of k for the case 
K = 0. The bottom figure is the top view of the upper one. For this singular 
case we notice the two lines at fc = 1/2 corresponding one to collapse and one 
to infinite expansion. On the projection plot (bottom), the trajectories k < 1/2 
and A; > 1/2 seem to intersect on those two fines. For A > cofiapse is 
not possible, and the trajectories do not intersect. We notice that close to the 
collapse condition k = 1/2, the values k < 1/2 describe a collapse motion, while 
fc > 1/2 describes the expansion motion. 
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Figure 13: The effective potential V(X) in the singular case fc = 1/2 = 0) 
for the critical energy corresponding to A = 1/4: the double root case. 
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Figure 14: Evolution of the effective potential V{X) in the singular case k = 1/2 
{5 = 0) within the energy range corresponding to possible collapse A e [1/2, 1]. 
We notice that unbounded motion occurs. 
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Figure 15: Period of the periodic motion (A < 1/2) motion versus |A — 1/2| 
in the singular case k = 1/2 {6 — 0). We notice a power-law divergence of the 
period as we approach collapse. The measured exponent is 3/2. See appendix 
^ for the analytical computation. 
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Figure 16: Phase space portrait for different values of A close to Ac = 1/2, in 
the case K = 1 and A; = 1/2 (5 = 0). The bottom figure is the top view of the 
upper one. We notice the singular evolution of the motion to larger scales as A 
approaches Ac. 
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Figure 17: The approach to cohapse. Initial positions available for the vortex 
with negative strength. The two vortex with strength 1 are fixed, and located by 
the two thick points on the plot. Circles refer to different values of the constant 
K, while the oval curves correspond to different values of A. The shaded region 
corresponds to initial conditions leading to the aperiodic infinite expansion type 
of motion. 
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Figure 18: Vortex trajectories during collapse. Initial positions are shown by 
circles, they correspond to a fastest collapse value A = a/3/2; negative vortex 
trajectory is plotted in black, positive ones in gray. All three vortices converge 
to the center of vorticity in finite time r = Att/3. 
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Figure 19: Unbounded near-collapse motion for k = 1/2. Negative vortex: 
black; positive: gray. Closest approach positions are marked by asterisks. 
Upper plot K < 0: vortices contract until a collinear configuration with Xmin = 
K/{l/2-A) is reached. 

Lower plot K > 0: closest approach is this case is Xmin ~ ~ A) and 

corresponds to an isosceles triangle. 
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Figure 20: Motion in the vicinity of a saddle-point A = Kc^ Contracting vortices 
approach an unstable equilibrium configuration, their spiral trajectories tighten 
and become nearly circular. 

Upper plot A < Kc^'- vortex configuration cannot pass through the saddle-point, 
and reflects back after reaching collinear conflguration. 

Lower plot A > Kc^: vortices "above" the saddle-point, and continue contraction 
(only a small piece of this part of the trajectory is shown). 
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Figure 21: Far from collapse S = 0.4. Poincare section of 253 passive tracers in 
the flow field generated by three point vortex. The run is over 4000 periods. The 
constant of motions are A = 0.9, K = 0. The vortex strengths are (—0.1, 1, 1). 
The period of the motion is T = 6.59. 
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Figure 22: Far from collapse S = 0.09. Poincare section of 253 passive tracers in 
the flow field generated by three point vortex. The run is over 4000 periods. The 
constant of motions are A = 0.9, K = 0. The vortex strengths are (—0.41, 1, 1). 
The period of the motion is T = 36.86. 
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